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Introduction 


The phenomenon of peeling and debonding of thin layers is a subject of interest 
to those concerned with adhesives, thin films, and layered materials. In recent years 
much attention has been focused on such problems as a result of increased interest and 
application of advanced composites and thin film coatings. (See for example ref. 1. 
An extensive list of references pertaining to the subject can be found therein ) A 
related problem which is of interest for its own sake but also represents a simple 
example of a tangled adhesive strip and of coplanar delamination interaction is the 
problem of a looped adhesive strip. This is the subject of the present study. 

We consider here the problem of an elastic strip which possesses an adherend on 
(at least) one of its surfaces. If the strip is deformed so that two portions of 
such a surface are brought into contact, a portion of the strip becomes bonded and a 
loop is formed (Figure 1). We shall be interested in determining the equilibrium 
configuration of such a strip and investigating the behavior of the strip when its 
edges are pulled apart. 

The problem shall be approached as a moving interior boundary problem in the 
calculus of variations with the strip modeled as an inextensible elastica and the 
bond strength characterized by its surface energy.* A Griffith type energy 
criterion shall be employed for debonding, and solutions corresponding to the problem 
of interest obtained. The solution obtained will be seen to predict the interesting 
phenomenon of "bond point propagation", as well as the more standard peeling type 
behavior. Numerical results demonstrating the phenomena of interest are presented as 
well and will be seen to reveal both stable and unstable propagation of the boundaries 
of the bonded portion of the strip, depending upon the loading conditions 
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*Bottega, W.J.: Peeling and Bond Point Propagation in a Self-Adhered Elastica 

To appear in Quart. J. Mech. and Appl. Math. 
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Formulation of the Problem 

Consider a thin elastic strip which possesses an adherend on one of its surfaces, 
and let the strip be closed on itself in a symmetric manner such that there exists a 
region in the plane of symmetry where the strip is bonded to itself (Figure 1). 
Additionally, let the edges of the strip be subjected to equal and opposite forces as 
shown. As a result of the inherent symmetry of the problem, only half of the strip 
need be considered in the ensuing analysis (Figure 2). The strip thus consists of a 
lifted segment, a bonded segment, and a looped segment. In what follows, all length 
scales have been normalized with respect to the half length, L, of the entire strip. 

We shall identify each point on the centerline of the strip by its normalized arc 
length, s, measured from the edge at which the external force is applied. In so doing, 
the half strip will be divided into four regions; corresponding to the lifted segment 
defined on 0 < s < a, a < s < b corresponding to the bonded segment , with the looped 
segment of the strip divided into two regions, defined by b < s < s* preceding the 
associated inflection point, s*. and s* < s < 1 following the inflection point. We 
shall be interested in assessing the behavior of the above system as a function of the 
magnitude of the applied load or the corresponding edge displacement. 

Let us first define the right handed cartesian coordinate system (x,y) as shown 
in the figure. In addition, let us define the angle 0(s) which measures the angle that 
the tangent of the strip at point "s" makes with the x-axis as s increases (see Figure 
2). One then may easily find the relations 



Figure 2 
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Energy Functional 


As we shall approach the problem as a moving interior boundary value problem in 
the calculus of variations, we next define the energy functional FI given in Table 1 
where U corresponds to the normalized strain energy of the segment of the strip 
defined on the domain D i , and is seen to be comprised solely of bending energy as 
the elastica is assumed to be inextensible . The strain energy of the (perfectly) 
bonded segment of the strip, i.e., the portion of the strip on [a,b] , is thus seen 
to vanish identically. The functional W corresponds to the normalized work done by 
the applied load. In that expression, P = PL 2 /D corresponds to the normalized 
counterpart of the magnitude of the applied load, P, and D represents the bending 
stiffness of the strip. The functional T corresponds to the "delamination energy" 
where 7 - yL /D represents the normalized counterpart of the surface energy of the 
bond, 7, while a 0 and b 0 correspond to the initial values of a and b respectively 
Finally, we introduce the constraint functional A, with Lagrange multiplier A 
which constrains the deflections of the segments of the strip on D. and D, to be 


continuous at s 
9 by eq. (i-a) . 


s*. The functions x 3 (s*) and x 4 (s*) may be expressed in terms of 
Thus , 


A " * / b cos ^ ds + A J ^ cos B ds . (2e') 

We note that the inclusion of A is equivalent to treating the segments on D 3 and D, 
as separate structures and including the work done by the internal force A , at s - s* 
It may be seen that the line of action of this force must be parallel to the x-axis 
as a result of the support condition at s - 1 (see Figure 2 ). 


Table 1 


4 r\ 

n = z LT - W + r - a 
i=l 

(1) 

where: 


U (,) =J |(d«/ds) 2 ds 

D. 

1 

(2a) 

0 

hi 

(2b) 

a 

W = -P x (o) = P / cos e ds 
o 

(2c) 

r = 2 7 (a-a o ) - 2 7 (b-b o ) 

(2d) 

A = a[x 3 (s*) -X 4 (s*)] 

(2e) 
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The Governing Equations 


The governing differential equations, boundary conditions, matching conditions and 
transversality conditions for the problem of interest are found by invoking the 
principle of stationary potential energy as shown in Table 2 below. In equation (3), 
the parameter 5 represents the variational operator. 

The transversality conditions (7a,b,c) result from the moving boundaries during 
peeling or bonding of the strip and thus are associated with equilibrium configurations 
of the system during these processes. 

The intermediate boundaries a, b and s*, as well as the inflection point angles 
a and a* are found as part of the solution to the problem. 

Table 2 


Principle of Stationary P.E.: 5(1 - 0 ^ 

(1) into (3): 


hT^) 2 + P-(cos a cos 9 .) = 0 , (i = 1,3,4) (4) 

2 v ds' i i i 


e 2 = */2 


where: i.= is*-* , “ p 05 e (0) , P P • 

(5) 


"2,3,4 9 ’ a 3,4 

a* = e(s*) . 

' P 3,4* A 


with B.C.s and M,Cs: 


and T.C.s: 


^(a) = * 3 (b) = w/2 , 

(6a, b) 


(7a) 

* 4 0) = 0 

(6c) 

G < b ‘>' ^ 2| s=b* = 2y 


* 3 (s*) = « 4 (s*) = a* 

(6d) 

(7b) 

1 

jf costf ds = -J *cos e ds 

D ^ 

(6e) 

ds, d« 4 

-^1 =-* | * = 0 

ds 's = s* ds s = s* 

(7c) 
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Criteria for Propagation of the "Bond Zone" Boundaries 


The conditions (7a) and (7b) establish the bond zone boundaries during their 
propagation and state that the values of a and b corresponding to equilibrium config- 
urations of the system during propagation of each interior boundary are those for 
which the bending energy densities at the point s - a' and s - b + are just balanced 
by the energy of the bond. In this context the quantities G(a') and G{b + } may be 
identified as the 'energy release rates" at the bond zone boundaries. The above 
suggests the criteria for propagation of the boundaries of the bonded region of the 
elastica, as listed in Tables 3a and 3b. 

Peeling 

If, for some initial a = a Q , eqn. (8a) is satisfied, no peeling will occur and a 
will remain at its initial value a Q with the lifted segment bending away from the 
plane of symmetry. If, for some initial a = a OI eqn. (8b) is satisfied, the lifted 
segment of the strip peels away from the plane of symmetry such that the value of a 
increases until the corresponding equality (7a) is satisfied. 

Following the above reasoning, we conclude that for the loop to maintain its 
initial configuration, conditions must be such that eqn. (9a) is satisfied. If 
eqn. (9b) holds, the loop would open as a result of excess bending energy at its 
e( 3& e > with b taking on smaller values until the energy of deformation is just 
balanced by the energy of the bond. 


Table 3a 


Criteria for Propagation 
of Bond Zone Boundaries: 


Peeling: 

If G{a„) < 2-i 

No Peeling 

(8a) 

If G{a J > 2 t 

Peeling until a 
satisfies equality (7a) 

(8b) 

Similarly, 

If G{b„'} < 2-i 

Loop maintains initial 
configuration 

(9a) 

If G{b.*> > 2y 

Loop grows (b decreases 



until b satisfies equality (7b) 

(9b) 
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Bond Point Propagation 


If conditions are such that the bending energy G{a } is sufficiently large while 
G(b 0 + } is sufficiently small, the bond zone boundary "a" will increase while b remains 
at its initial value until a - b". At this stage, as the resultant bending moment 
m ba acting at the bond point satisfies eqn (10) (Table-3b), and hence acts in a 
clockwise sense, tending to rotate the strip in this sense, while simultaneously 
there exists a sufficient surplus of bond energy to counter the bending energy of the 
loop and induce bonding at the loop edge, the strip behaves locally as if rolling 
over its counterpart and the ’’bond point” s * a - b will propagate such that the loop 
closes and shrinks until the corresponding equality (7b) is satisfied. At this point 
the surplus bond energy is depleted and bonding at s = b + can no longer occur. 
Equivalently, the growth condition (9b) as well as (8b) will become satisfied and 
conditions will then be such that propagation can occur in both directions simul- 
taneously. Under such conditions, the loop will expire and the strip will separate 
completely . 


Table 3b 

"ROND POINT PROPAGATION”: 


If (9a) and (8b) satisfied, a -*• b o . 


Resultant Bending Moment: 


m ba = [[d»/ds]] s = b = [d«/ds] s=b * - [d«/ds] s=b - < 0 


(thus, 


m^is clockwise) 


( 10 ) 


Thus, "bond point" propagates until equality (7b) satisfied 
and loop expires. 
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Basic Integrals 


The nonlinear differential equations (4) governing the local rotations on each 
segment of the strip are seen to be of standard form, which upon solving for ds 
integrating over the corresponding domain, and using the transformation given by ’(11) 
results in expressions for the segment arc lengths (i-1,3,4) in terms of the 
inflection point angles and corresponding external or internal loads as given in 
Table 4. We note that these lengths will vary as a result of bonding and debondine. 
The functionals F(q,$) and F k (q) in eqns. (12) correspond to elliptic integrals of 
the first kind and complete elliptic integrals of the first kind, respectively 
defined by J ' 

$ d <t>. 

F( V $) -/o ■ and F k ( q i ) " F (<*i.*/ 2 > • (ii-a,b) 

yi-qj^ sin ^ 

We shall first consider equilibrium configurations of the lifted segment 

(fleDj) and of the looped segment (0eD 3 +D<) separately, and then examine their 
interaction. 


Table 4 

BASIC INTEGRA! ,S 
Transformation: 

sin(«. /2) = q. sin </> . (11a) 

q; = sin («./2) (i = 1,3,4) (lib) 

(11) into (4), solving for ds and integrating =» 
segment arc lengths : 

1 1= a = [F^qj) - F(q 1 >* 1 )/#’ (12a) 

« 3 = S* - b = [F k (q*) - F(q*,#3>]/yr (12b) 

f 4 = 1-s* = F k (q‘)/yi (12c) 

where q 3 =q 4 =q* (13a) 

* j = sin ' 1 {l/(q i ^>> (13b) 

F(q,§) - Elliptic Integral of 1st Kind 
Fk(qj) - Complete E.I. of 1st Kind 
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Behavior of the Lifted Region 


The deflection of the edge of the strip at which the load is applied is found 
by solving equation ( 4 a) for d 0 x /ds and substituting the resulting expression into 
equation (i-a) , with appropriate limits of integration. Then, upon incorporating 
the transformation (13), we obtain the load edge deflection, A 0 , as given in Table 5, 
where 

E(q i ,<l>) - yi-q i 2 sin 2 ^ i d^ and E^q^ ” E(q i> 7r/ 2 ) » (iii-a,b) 

correspond to elliptic integrals of the second kind, and complete elliptic integrals 
of the second kind, respectively. The explicit form of the transversality condition 
at the "trailing edge" of the bonded region (i.e., at s - a) may be found by solving 
equation ( 4 a) for [d^ 1 /ds] s; , a and substituting the resulting expression into ( 7 a). 

We then have the condition which (implicitly) defines the location of the trailing 
edge of the "bond zone", during peeling, given by eqn. (15). 

Substitution of equation (15) into equation (12), with i = 1, and equation (14) 
gives explicit relations for the magnitude of the applied load as a function of "a" 
and the normalized load point deflection as a function of "a" respectively, with the 
load point rotation a a parameter. Specific results corresponding to selected values 
of 7 will be presented in a later section. 

Table 5 

BEHAVIOR OF LIFTED REGION 
LOAD POINT DEFLECTION: 

A o S ' x| s = o ~ “« 2E <VV ' F ( q l’ # p] 


- [2E k ( qi ) - F k ( qi )]} (14) 


where E(q,*) - Elliptical Integral of 2nd Kind 


E. (q .) - Complete E.I. of 2nd Kind 


T.C. @ s = a’ (7a) becomes 

P = - 27 /cos a (jt/ 2 < a < jr) (15) 
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Equilibrium Configurations of the Loop 


The angle which measures the rotation of the tangent of the loop at its 
inflection point s - s* is found by imposing the condition (6e). Thus, solving 
equations (4b) and (4c) for dfl/ds, substituting the resulting expressions into the 
left and right hand sides of (6e) respectively, noting (6d) and incorporating the 
transformation (11) for i - 3,4 we obtain the condition given by (16) below where 
as defined earlier, a* - c* 3 - a 4 and q* — q 3 - q*. It may be seen from equation 
(16) that a* is independent of the size of the loop, of the material and geometric 
properties of the strip and bond, and of the magnitude of the applied load, and 
thus is a "characteristic angle" of the problem. Equation (16) may be solved 
numerically to yield the value of a* as given below. 

The total (half) arc length of the loop, i, is simply comprised of the sum 
of the lengths of its constituent segments. Thus, adding equations (12b) and 
(12c) yields the relation for £, given by (18). The relative portions of the loop 
corresponding to its constituent segments are then found by dividing eqns.(12b) 
and (12c) by (18). We thus have 

£ 3/ £ ~ [ F k (q*) - F(q*,$*) ]/[2F k (q*) - F(q*,$*)] (19a') 

and 

i-t/z - F k (q*)/[2F k (q*) - F(q*,$*)] , (19b') 

which are seen to correspond to "characteristic length ratios" of the problem. The 
above ratios may be evaluated, using the computed value of a*, to yield the values 
given at the bottom of Table 6. 


Table 6 

EQUILIBRIUM CONFIGURATIONS OF THF. LOOP 
Imposing (6e): 

2[2E k (q-> - F^q*)] = 2E(q*, .*) - F(q*, #•) (16) 


(where #* = *^ = # ). 

Solving (16) yields "characteristic angle" of inflection point 

a* = 117.54* (for any loop size and 

mat’l./geom. props.) (17) 

(12b) + (12c) =» loop (half) arc length t : 

1 = * 3 + < 4 = [2F^q*) - F(q*. (18) 

(12b,c)/(18) =» "Characteristic Arc Length Ratios": 
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tj/« = 0.3254 and t./t = 0.6746 

(for any loop size and mat’l./geom. props.) 


(19a,b) 


Debonding of the Looped Segment 


We may next consider equilibrium configurations of the looped segment of the strip 
as it opens (debonds) by evaluating the explicit form of the transversality condition 
at s * b + in a manner analogous to that done earlier for the lifted segment. Doing so 
we find that during opening of the loop, the condition ( 7 b) takes the form given by 
eq'n.(20) where a* is given by (17). Since, as discussed earlier, a* is a "character- 
istic angle" (i.e. , it maintains a fixed value for any equilibrium configuration of the 
loop,) it is seen from the above expression that the internal force A maintains a 
constant value during debonding of the loop. 

Substitution of (20) into equation (18) yields a critical value of the loop length 
given by equation ( 21 ) below, with the inequalities ( 9 a, b) now interpreted in terms of 
the arc length of the loop; e.g. - ( 22 a, b). 

It may be noted that a minimum value of the normalized bond energy is required for 
the elastica to remain adhered to itself. This value 7 ■* 7 min corresponds to the 
limiting case where the loop traverses the entire strip and is in self-contact only at 
the loading edge s = 0 (i.e., it corresponds to the limiting case when JL cr — 1). Upon 
employing the result (19) we find the desired value given below. Adherends whose 
normalized bond energies possess magnitudes which are below this value are thus not 
"strong enough" to maintain a self-adhered configuration. 

Table 7 

T.C. @ s = b + (7b) becomes: 


a = - 27 /cos a* (= constant for given 7) (20) 


(20) into (18): 


< = [2F j^q*) - F(q*, #*)]y(-cosa*)/2 7 (21) 


If t > t 

cr 

If t < t 

cr 


No debonding of loop occurs 

Debonding of loop occurs 
1 increases (b decreases) until t 


t 

cr 


= 1-7 


7 min 


= 2.292 


(22a) 


t 

cr 


(22b) 


(23) 
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Effective Bond Strength and Propagation Behavior 


A plot of Z cr versus 7 is displayed in Figure 3. It may be observed from the 
figure that the amount of "effective bond strength" gained, as measured by the relative 
decrease in i cr> significantly decreases as 7 is increased beyond 200 . 

Peeling and Bond Point Propagation 

The above solution offers the following scenario for a looped adhesive strip with 
given 7 > 7 min existing in an initial configuration such that Z Q - l-b Q > Z cr (or 
equivalently b G < b cr « l-i cr ), with an initial lift zone size of - a D . As P is 
increased the corresponding value of — a is increased according to equation ( 12 a ) , 
with the associated deflection A c varied according to equation (14). This process is 
continued until equation (15) is satisfied at which point peeling begins with the 
"trailing edge" of the "bond zone" s = a propagating (and Z x increasing). As the 
initial loop length Z Q is larger than the critical length, the loop edge boundary of 
the bond zone remains at its initial value until a - b c . At this point if conditions 
are such (equation (15) satisfied) that peeling continues, the bond point a = b 
propagates with Z decreasing (b increasing). During this phase the values of a*, 

^ 3 /^ » and Z^/Z maintain the values given by (17) and (19a,b) respectively. The loop 
thus shrinks in size during this phase, with its geometry evolving through successive 
self-similar shapes (as if the strip were being pulled through a rigid clamp at s - b). 
This process continues until Z - Z cx , at which point conditions are such that peeling 
may occur at both s ■* b and s “ b (i.e, , in both directions) simultaneously. At 
this instant the loop is terminated and the surfaces on each side of the plane of 
symmetry separate. Results corresponding to specific values of 7 are presented 
next. 
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Fig. 3. Variation of Critical Loop (Half) Length with Normalized Bond Energy 
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Numerical Results 


Results are presented corresponding to selected values of the normalized bond 
energy. Specifically, we consider the strip/adherend systems whose material and 
geometric properties are characterized by the values 7 — 5, 10, 50, 100, and 1000 
respectively. 

The projections of the associated lift zone/bond point "propagation paths" 
in the A 0 -a, P-a and P-A 0 spaces are calculated and are displayed in Figures 4-7. 

Each curve is terminated at the critical values a - b cr = l--* C r> which are given 
by the values b cr - .3230, .5213, .7859, .8486, and .9522 for the respective values 
of 7 considered. The prepropagation load-deflection behavior of lift zone segments 
corresponding to initial lengths of a 0 = 0.25 and a Q - 0.50 are also displayed in 
Figure 6 . Finally, the variation of the magnitude of the internal force A , as a 
function of the loop length 2 , is displayed in Figure 8 . 

It may be seen from the figures that propagation of the lift zone boundary or 
bond point occurs in a stable manner for a deflection controlled test, and in an 
unstable and "catastrophic" manner for a force controlled test. The following example 
illustrates the general behavior of the self-adhered elastica. 



Fig. 4. Lift Zone/Bond Point Propagation Paths (7 - 5, 10, 50, 100, 1000): 
Load Edge Deflection vs. Bond Zone Boundary 
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Let us consider a strip/adherend combination with 7 - 50 which is initially 
configured such that a 0 - 0.25 and b 0 - 0.60. It may be observed that a system 
characterized by one of the smaller values of 7 considered could not maintain a loop 
size this small (b 0 this large) and that for such a system the loop edge boundary b 
would immediately decrease to its corresponding critical value. Returning to the’ ’ 
previously defined case (7 - 50) , the system initially follows the prepropagation 
path for a Q - 0.25 in Figure 6 , as P and A 0 increase from the origin. During this 
phase, the system simultaneously follows a purely vertical path, corresponding to 

3 " ? - ’ 1 " . the A °' a and P ' a s P aces (Figures 4 and 5). The system continues to behave 

m this fashion until the propagation path corresponding to 7 - 50 is intercepted. 

At this point the lifted segment of the strip has accumulated enough bending energy at 
the bond zone boundary s - a , for the lift zone to propagate. We note that as 
o > ^cr < b cr ) the bending energy of the looped segment at s = b + is insufficient 
for propagation so b remains at its initial value. Let us first consider the case 
where the load edge deflection, A OI is controlled. For this case, as A is 
incrementally increased, a increases incrementally following the corresponding path 
in Figure 4. The corresponding values of P may be observed, from Figs. 5 and 6 to 
decrease accordingly. This process continues, with the strip peeling from its 
symmetric counterpart in a stable manner, until a - b D . At this point the bending 
energy of the loop at s - b is still sufficiently low as to maintain the bond, 
while that at s - b is large enough to permit debonding. As A 0 is increased 
further, the bond point s - a - b then propagates in a stable manner, with the loop 
shrinking through a series of self- similar shapes until Z = Z cz , at which point 
sufficient bending energy exists on both sides of the bond point and the strip 
separates. For the case where P rather than A 0 is controlled, the system behaves 
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Fig. 5. Lift Zone/Bond Point Propagation Paths (7 — 5, 10, 50, 100) • 
Applied Load vs. Bond Zone Boundary 
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in an analogous manner except that all behavior during the propagation phases occur in 
an unstable manner. Thus, for this case the process of complete separation of the 
elastica is initiated as soon as P reaches a critical value. 

The phenomena described above may be observed by simply closing a piece of 
adhesive tape on itself, thus forming a loop, and then peeling the edges apart. Such 
an "experiment” would correspond to a deflection controlled test, with the normalized 
bond energy characterizing the tape observed to be at the upper end of the range 
of values considered in the numerical simulations presented herein. 



Fig. 6 . Prepropagation and Lift Zone/Bond Point Propagation Paths (7 - 5, 10, 
50, 100): Applied Load vs. Load Edge Deflection 


it) 00 



Fig. 7. Lift-Zone/Bond Pt. Propagation Paths (7 = 1000): (a) Applied Load vs. Bond 

Zone Boundary, (b) Applied Load vs. Load Edge Deflection 



Fig. 8 . 



Internal Force at Inflection Point of Loop as a Function of 

Loop (Half) Length 
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